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Abstract 

Non- perturbative exact flow equations describe the scale dependence of the 
effective average action. We present a numerical solution for an approximate form 
of the flow equation for the potential in a three-dimensional A-component scalar 
field theory. The critical behaviour, with associated critical exponents, can be 
inferred with good accuracy. 


Exact non-perturbative renormalization group equations describing the scale depen¬ 
dence of some type of effective action have been known for a long time [|T| . They account 
for the consecutive inclusion of fluctuations in a held theory. There exist many versions 
of such equations which are, if exact, all equivalent, since they all describe in one way or 
another properties of the (Euclidean) functional integral which defines the theory. The 
most difficult part, however, is not so much to derive an exact how equation - this usually 
follows from simple manipulations of the functional integral. The challenge is rather to 
find a formulation which can be used for practical computations beyond perturbation 
theory. For such practical purposes the most general form of the effective action has 
always to be truncated. Solutions of the truncated equations are only approximations to 
the exact how equations. If a small parameter is available one can often organize a series 
of truncations as a power series in the small parameter - this is how usual perturbation 
theory is recovered. If no small parameter is known a priori one has to use as much as 
possible the knowledge of properties of the model in order to conceive a useful trunca¬ 
tion. This is where it becomes important what effective action is selected and what the 
particular form of the how equation is. Only if the effective action has a simple physi¬ 
cal meaning so that its couplings can also be understood by other methods, and if the 
evolution equation has a form which incorporates directly the most prominent properties 
of the system, there is a realistic chance of getting a working non-perturbative tool. Be¬ 
yond perturbation theory exact how equations should, therefore, not be viewed as mainly 
a consistent (but unsolvable) mathematical system, but rather as a convenient starting 
point for appropriate non-perturbative expansion methods. 

Typically, the truncated how equations constitute a system of non-linear partial dif¬ 
ferential equations which can be solved by numerical methods. It is the aim of this letter 
to report on the development of algorithms which are adapted to the specihc numerical 
problems of this approach. As an example we investigate here the IV-component Heisen¬ 
berg models in three dimensions. This is a classical testing ground for non-perturbative 
methods. We will extract the critical exponents of these models from a numerical solu¬ 
tion of the how equation for the scale dependent effective potential, as obtained from the 
lowest order in a derivative expansion of the effective action. 

Recently the average action T^ [[2| has been proposed as the free energy with an infrared 
cutoff ~ k. It is formulated in continuous space so that all symmetries of the model are 
preserved. More precisely the average action includes the effects of all huctuations with 
momenta q 1 2 > k 2 . In the limit k —► 0 it becomes the standard effective action (the 
generating functional of the 1PI Green functions), while for k —» oo it equals the classical 
action. There is a simple functional integral representation [2] of T k also for k > 0 
such that its couplings can, in principle, also be estimated by alternative methods. The 
average action is a coarse-grained free energy in the sense that short distance fluctuations 
are already integrated out. [] The exact non-perturbative flow equation for T^ takes the 

1 One should not use IT directly to replace the action in the functional integral, as this is the role of 

the “cutoff action” used in earlier approaches jT]. For this purpose an explicit ultraviolet cutoff term has 
to be added in order to prevent double counting of the short distance modes. This cutoff term is easily 
computed in simple theories of scalars and fermions [^) , much more involved for an abelian gauge theory 
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simple form of a renormalization group improved one-loop equation [[| 

d t r k = k^-T k = -Trd t In (rf + R k ) , 

where d t acts only on the infrared cutoff piece R k as 

d t In (r{ 2) + R k ) = (r< 2) + R fc ) _1 k^R k . (2) 

The trace involves a momentum integration and summation over internal indices. Most 
importantly, the relevant infrared properties appear directly in the form of the exact 
inverse average propagator r(, , which is the matrix of second functional derivatives with 
respect to the fields. There is always only one momentum integration - multi-loops are 
not needed - which is, for suitable R k , both infrared and ultraviolet finite. 

The flow equation ([TJ) seems to be a good starting point for a non-perturbative ap¬ 
proach. Nevertheless, it remains a complicated functional differential equation without 
any chance to be solved exactly. Approximate solutions need truncations - and the crux of 
the problem lies there. One possibility is to keep only a few invariants in T k and, thereby, 
reduce eq. (jl|) to a finite set of ordinary differential equations for a finite number of 
couplings. Very satisfactory results have already been obtained this way ||. Much more 
information is contained if, instead of a finite number of couplings, arbitrary functions 
of one or several parameters are considered. Examples are an arbitrary field dependence 
of the average potential U k (p) in scalar theories (with p = \<$> 2 ) § 1, or an arbitrary 
momentum dependence of the two-point function G(p ) []7| or even the four-point function 
A(gi, ^ 2 , <? 3 , Qa) !• The functional differential equation becomes then a partial differential 
equation for a function of at least two variables, e.g. U(p,t ) = U k (p), or a system of 
partial differential equations for several such functions. If one does not want to resort 
to further approximations at this level, one needs appropriate tools for solving this type 
of partial differential equations. Analytical solutions can be found only in certain lim¬ 
iting cases and for most purposes numerical solutions seem the adequate tool. This is 
not so straightforward as it may seem at first sight, since the solutions of the differential 
equations have in general a highly unstable character due to the presence of relevant pa¬ 
rameters. Fine tuning of initial conditions is needed in order to be near a phase transition 
- or, in particle physics language, to have renormalized masses much smaller than the 
ultraviolet cutoff. The algorithm for a numerical solution must be compatible with this 
situation and guarantee numerical stability of the critical solution which corresponds to 
the phase transition. 

In this letter we demonstrate the capacities of such a non-perturbative method by 
computing the fc-dependent average potential U k (p) for an A r -component scalar field the¬ 
ory directly in three dimensions. For k —> 0 this gives the free energy which encodes the 
equation of state of the system. In the phase with spontaneous symmetry breaking the 
minimum of the potential occurs for k = 0 at p 0 7^ 0- The massless Goldstone excita¬ 
tions around this minimum are notoriously difficult to treat by alternative methods. For 

H , and not yet found for non-abelian gauge theories. 


2 



example, the standard loop expansion is order by order highly infrared divergent. In the 
symmetric phase the minimum of U k {p) ends at p 0 = 0 for k — 0. The two phases are 
separated by a scaling solution for which Uk/k 3 becomes independent of k once expressed 
in terms of a suitably rescaled field variable p. 

Our truncation is the lowest order in a systematic derivative expansion of i, |,i 

T fc = J d d x{u k (p) + ^Zkd^ad.r}- (3) 

Here 0“ denotes the TV-component real scalar field and p = ~0 a 0 a . We keep for the 
potential term the most general O(AQ-symmetric form U k (p), whereas the wave func¬ 
tion renormalization is approximated by one ^-dependent parameter. Next order in the 
derivative expansion would be the generalization to a p-dependent wavefunction renormal¬ 
ization Zk(p) plus a function Y k (p) accounting for a possible different index structure of 
the kinetic term for N >2 i,0- Going further would require the consideration of terms 
with four derivatives and so on. For the three-dimensional scalar theory the anomalous 
dimension p is known to be small and the derivative expansion is, therefore, expected to 
give a reliable approximation [5], The main reason is that for rj — 0 the kinetic term in 
the ^-dependent inverse propagator must be exactly proportional to q 2 both for q 2 —> 0 
and q 2 —> oo. This holds for arbitrary constant “background” field (j) a . Similar, although 
less stringent, arguments concern the smallness of the p-dependence of the kinetic term 
||. For the scaling solution for N — 1 this weak p-dependence has been established 
explicitly ||. We finally mention that p is proportional to a small parameter A /87 r 2 , 
where A is a suitably defined quartic scalar coupling. (For the scaling solution A takes 
a fixed point value A*/ 87 r 2 = 0.12 for N = 1.) We expect that the derivative expansion 
can be understood as an expansion in this small parameter. However, since some of the 
other parameters effectively behave ~ A -1 this expansion is not equivalent to the usual 
perturbative polynomial series in A. 

For a study of the behaviour in the vicinity of the phase transition it is convenient to 
work with dimensionless renormalized fields f] 

p =Z k k 2 ~ d p 

u k (p) =k~ d U k {p). (4) 

With the truncation of eq. (||) the exact evolution equation for u' k = du k /dp [||] reduces 
then to the partial differential equation 

9" = (-2 + vK + (d- 2 + v)K 

- 2v d (N - l)u k lf(u' k ; p) - 2v d (3u k + 2 pu k )lf(u' k + 2pu k ; p), (5) 

where t = In {k/ A), with A the ultraviolet cutoff of the theory. The anomalous dimension 
p is defined by 

d = ~^ t lnZk ( 6 ) 

2 We keep the number of dimensions d, arbitrary and specialize only later to d = 3. 
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and 


v,I 1 = 2 d+ Vir , (7) 

with t ’ 3 = 1 / 87 T 2 . The “threshold” functions l d (w,r)) result from the momentum integra¬ 
tion on the r.li.s. of eq. and read for n > 1 , with y = q 2 /k 2 

i*(w;rj)= ~ n j 0 d yy^ + 1 -^-[y ( 1 + r (y)) + w ]~ (n+1) 

- 7:V [ dyyir(y) [y( 1 + r(y)) + w]~ (n+1) . (8) 

2 jo 

Here r(y) depends on the choice of the momentum dependence of the infrared cutoff and 
we employ 

r(y) = -■ (9) 

1 — e v 

This choice has the property lim g 2^ 0 R k /Z k k 2 = lim y _, 0 yr(y) = 1, whereas for q 2 k 2 
the effect of the infrared cutoff is exponentially suppressed. The “threshold” functions 
account for the decoupling of modes with mass larger than k and decrease rapidly for 
w 1. For our purpose we use numerical fits for these functions. Finally, the anomalous 
dimension is given in our truncation by 

V(k) = ^^K\ 2 m d 2a (2\K), ( 10 ) 

with k the location of the minimum of the potential and A the quartic coupling 


u' k {n) =0 

u'I(k) =X. (11) 


The function m\ 2 is given by 


m 2 } 2 (w ) = / dyyi 


dr 


2-2 


1 + r + y 


dr 

dy 


(1 + r ) 2 [(1 + r)y + w\ 
2 / 


d 


dr\ dr 


2?/7r + 2 h/^ r-2y 1 + r + y— — 


dy 


dy 


dy I dy 


+ 


(1 + r)y (1 + r)y + w 


( 12 ) 


We point out that the argument 2\k turns out generically to be of order one for the 
scaling solution. Therefore, k ~ A -1 and the mass effects are important, in contrast to 
perturbation theory where they are treated as small quantities ~ A. 

Our aim is the development of algorithms for the numerical solution of the partial 
differential equation (|5|), and a comparison with previously used expansion methods. Near 
the phase transition the trajectory spends most of the “time” t in the vicinity of the critical 
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k- independent scaling solution given by d t u'*(p) = 0 . [] Only at the end of the running 
the “near-critical” trajectories deviate from the scaling solution. For k —> 0 they either 
end up in the symmetric phase with k = 0 and positive constant mass term m 2 such 
that i4(0) ~ m 2 /k 2 ; or they lead to a non-vanishing constant po indicating spontaneous 
symmetry breaking with k —> Z 0 k 2 ~ d p 0 . The equation of state involves the potential 
Uo(p) for temperatures away from the critical temperature. Its computation requires the 
solution for the running away from the critical trajectory which involves the full partial 
differential equation (||). We have developed two alternative numerical approaches which 
we briefly describe in the following. For both methods we replace the variable p by a 
discrete set of points p;, i = 1, 

I) Consider first what happens if we make a Taylor expansion of Uk around some 
arbitrary point pi 

OO 1 

Mp) = J2 —u[ n \p - pi) n , ( 13 ) 

n =0 n - 

with u\ n \k) = u^\pi). The potential is then described by infinitely many couplings 
u[ n \k). The flow equations for these couplings are obtained from appropriate p-derivatives 
of eq. evaluated at p = p+ We observe that the flow equation for u\ V) involves u^\ 
uf 1 and uf\ the one for u ^ needs in addition uf\ and the system is never closed. 
Our first approach considers at every point pi the differential equations for u\ 1 ^ and uf\ 
The “missing” couplings and ] appearing in these equations are determined by 
matching the expansion around pi with similar expansions around different points pj^i- 
More precisely, the matching is done by equating fourth order polynomial expansions of 
u' k (p) around two neighbouring points p t and p i+ 1 at half-distance, and similarly for u'l(p) 


u 


k)i 


Pi + Pi +1 


_ (1) . (2) Pi+1 

= U) ’ + U) -- 


A + u f)A±L 


Pi 


+ u) 


(4) (,Pi+ 1 Pi) 
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Pi + Pi+l \ _ (1) (2) Pi+1 ~ Pi (3) {Pi+ 1 ~ Pi) 2 (4) {Pi +1 ~ Pi) 
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k)i 


Pi + Pi+1 
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+ U. 


( u k)i +1 


Pi + Pi+1 ^ _ (2) 


— U i +1 U i +1 


2 

(3) Pi+1 


i+1" 


Pi + ^(4) (Pi+1 


— U. 


i+1' 


48 


Pi 


| „,( 3 4 ) ( P ®+1 

"r U i +1 


Pi) 


(14) 


(15) 


Combining eqs. (|T4|) and (p~5D at all m — 1 intermediate points (i = 1, ...,m) gives an 
algebraic system of 2m — 2 equations for the 2m unknowns uf^ and uf 1 . In order to 
obtain the remaining two necessary equations, we also match for the initial and end 
points the third derivative of the expansions, at (pi+p 2 )/2 and (p m _i+p m )/2 respectively. 
Together these equations make up an algebraic system which has a unique solution. This 


3 The resulting ordinary differential equation for u*(p) has already been solved numerically for a 

somewhat different choice of the infrared cutoff 101. 
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allows one to express the couplings uf 1 and uP as functions of the couplings uP and 
Uj, j = 1 The integration of the remaining system of 2m differential equations 

for the couplings uf\ u\ 2 * is done with the fifth order Runge-Kutta algorithm using 
the embedded fourth order method for precision control. The function u' k (p ) is finally 
reconstructed by patching the fourth order Taylor expansions around pi together at half¬ 
distance around neighbouring points. The polynomial patching improves with decreasing 
distance between neighbouring expansion points. This can be used to check the stability 
of numerical results. 

We note that the matching conditions of eqs. (0), (ED guarantee continuity of u' k (p) 
and u'l(p). (Further smoothcning could be applied for up (p) if needed.) Furthermore, 
they imply that the differential equations for uP and uP do not only incorporate infor¬ 
mation from directly neighbouring points, but from the whole range of points, as implied 
by the algebraic solution for uf* and v!p. Nevertheless, one observes that the contribu¬ 
tions of up 2 '* to uf A) (j 7 ^ i) rapidly decrease with increasing |i — j\. The decoupling 
from distant points could be used to obtain an approximate expression for uf A * which 
becomes useful if a large number of expansion points is considered. For the computation 
of the critical exponents we expand here the potential around 10 points. As long as the 
minimum of the potential is located away from the origin we choose expansion points 
pi proportional to hi. If the minimum is at or very close to the origin we use instead 
expansion points corresponding to fixed values of p instead of fixed p. 

II) For the alternative approach we first consider a standard discretized version of eq. 


u. 


(i) 


(k + Ak) - uf* (k) 

A k 


= F 


uP(k),u\ z> (k),uP(k) 


.( 2 ) 


(3) I 


(16) 


The higher p-derivatives uf\k), u\ 3 \k) can be inferred from the differences of the values of 
u i±i{k) an d the integration seems straightforward. The main problem with this approach 
is the appearance of numerical instabilities, with the numerical solution becoming strongly 
oscillating after a few integration steps. These instabilites are suppressed if the r.li.s. of 
eq. (ED i s evaluated at k + Ak instead of k. This leads us to replace the partial differential 
equation (||) by the system of m algebraic equations 


uP(k + A k) — u A *(k) 
: A k 


= F 


u A *(k + Ak), uf } (k + Ak ), uf’(k + Ak) 


.( 2 ) 


(3b 


(17) 


with the higher p-derivatives expressed in terms of uP l (k + A k). [] At every step the 
m unknowns uP(k + A k) can be calculated with the Newton-Raphson method. Further 
improvement in the accuracy can be achieved by making use of the solution at k and 
k — A k for the calculation of the solution at k + A k according to 


3 uP (k + Ak)- 4 u\ L> (k) + u\ 1} (k-Ak) 

2Ak “ 


,(i)/ 




uP (k + A k ), uY’ (k + Ak ), uf 1 (k + A k) 


.( 2 ) 


(3b 


4 For the end points p\, pm, the higher derivatives can be obtained at the same level of accuracy in 
terms of w^ i3)4 and up m _ i jm _ 2 , m - 3 respectively. 
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( 18 ) 


The expression on the l.h.s. approximates the ^-derivative at k + A k with an accuracy of 
0(|A£;| 2 ). We use 60 points for the discretization of the variable p and a varying number 
of k steps (around 500) until stability of the results is obtained. 

The use of two algorithms for the integration of eq. (|j) provides a good check for 
possible systematic numerical uncertainties. The two methods give results which agree 
at the 0.3 % level and the difference is most likely due to the different use of fits for the 
“threshold” functions. We expect the numerical solution to be an approximation of the 
solution of the partial differential equation (^) with the same level of accuracy. This has 
to be compared with the uncertainty induced by the omission of the higher derivative 
terms in the average action. The latter is expected to be of the order of 77 , as we discussed 
earlier, and is the main source of error for the results presented in the following. 

In fig. 1 we present the results of the numerical integration of eq. for d = 3 and 
N — 1. The function u' k (p) is plotted for various values of t — ln(/c/A). The evolution 
starts at k — A (t — 0) where the average potential is equal to the classical potential 
(no effective integration of modes has been performed). We start with a quartic classical 
potential parametrized as 

«a(p) = a a(p- k a ). (19) 

We arbitrarily choose Aa = 0.1 and fine tune so that a scaling solution is approached 
at later stages of the evolution. There is a critical value K cr ~ 6.396 x ICE 2 for which the 
evolution leads to the scaling solution without ever deviating from it. For the results in 
fig. 1 a value slightly smaller than n cr is used. As k is lowered (and t turns negative), 
u k (p) deviates from its initial linear shape. Subsequently it evolves towards a form which 
is independent of k and corresponds to the scaling solution d t u'*(p) = 0. It spends a long 
“time” t - which can be rendered arbitrarily long through appropriate fine tuning of ka - 
in the vicinity of the scaling solution. During this “time”, the minimum of the potential 
u' k (p) takes a fixed value /c*, while the minimum of Uk(p) evolves towards zero according 
to 

Po(k) = kn*/Z k . (20) 

The longer u k (p) stays near the scaling solution, the smaller the resulting value of po(k) 
when the system deviates from it. As this value determines the mass scale for the renor¬ 
malized theory at k — 0, the scaling solution governs the behaviour of the system very 
close to the phase transition, where the characteristic mass scale goes to zero. Another 
important property of the “near-critical” trajectories, which spend a long “time” t near 
the scaling solution, is that they become insensitive to the details of the classical theory 
which determine the initial conditions for the evolution. After u k (p) has evolved away 
from its scaling form u'^(p), its shape is independent of the choice of Aa for the classical 
theory. This property gives rise to the universal critical behaviour near second order phase 
transitions. For the solution depicted in fig. 1 u k (p) evolves in such a way that its mini¬ 
mum runs to zero with u k ( 0) subsequently increasing. Eventually the theory settles down 
in the symmetric phase with a positive constant renormalized mass term m 2 = k 2 u k (0 ) 
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as k —> 0. Another possibility is that the system ends up in the phase with spontaneous 
symmetry breaking. In this case k grows in such a way that p 0 (k ) approaches a constant 
value for k —> 0 . 

The approach to the scaling solution and the deviation from it can also be seen in fig. 
2. The evolution of the running parameters 7 c(t), A (t) starts with their initial classical 
values, leads to fixed point values tc*, A* near the scaling solution, and finally ends up in 
the symmetric phase (k runs to zero). Similarly the anomalous dimension r](k), which is 
given by eq. ([ 10 ]), takes a fixed point value 77 * when the scaling solution is approached. 
During this part of the evolution the wave function renormalization is given by 


Z k ~ k~ v * (21) 

according to eq. (fi). When the parts of the evolution towards and away from the fixed 
point become negligible compared to the evolution near the fixed point - that is, very close 
to the phase transition - eq. (| 2 l|) becomes a very good approximation for sufficiently low 
k. This indicates that 77* can be identified with the critical exponent 77 . For the solution 
of fig. 2 (N — 1) we find tc* = 4.07 x 10~ 2 , A* = 9.04 and 77 * = 4.4 x 10~ 2 . 

As we have already mentioned the details of the renormalized theory in the vicinity of 
the phase transition are independent of the classical coupling Aa- Moreover, the critical 
theory can be parametrized in terms of critical exponents ||, an example of which is the 
anomalous dimension 77. These exponents are universal quantities which depend only on 
the dimensionality of the system and its internal symmetries. For our three-dimensional 
theory they depend only on the value of N and can be easily extracted from our results. 
We concentrate on the exponent u, which parametrizes the behaviour of the renormalized 
mass in the critical region. The other exponents are not independent quantities, but can 
be determined from 77 and v through universal scaling laws J 9 | . We define the exponent v 
through the renormalized mass term in the symmetric phase 

777 2 = — k 2 u' k ( 0) for k —> 0. (22) 

Z k dp 

The behaviour of m 2 in the critical region depends only on the distance from the phase 
transition, which can be expressed in terms of the difference of from the critical value 
K cr for which the renormalized theory has exactly m 2 = 0. The exponent v is determined 
from the relation 

m 2 ~ \5ka\ 2i/ = |aca — K C r\ 2u - (23) 

For a determination of v from our results we calculate m 2 for various values of tca near 
K cr . We subsequently plot ln(m 2 ) as a function of In | <5 av a I • This curve becomes linear for 
Sk a —> 0 and we obtain v from the constant slope. In the past the critical exponents of the 
0(W)-symmetric theory were calculated from truncated versions of the partial differential 
equation (||) ||. The strategy was to turn eq. © into an infinite system of ordinary 
differential equations for the coefficients of a Taylor expansion analogous to eq. ([13]) 
around the “running” minimum of the potential. This infinite system was approximately 
solved by neglecting p-derivatives of u k (p) higher than a given order. The apparent 




convergence of the procedure was checked by enlarging the level of truncation. We now 
have an alternative way of estimating the accuracy of this method. Our numerical solution 
of the partial differential equation ([5]) corresponds to an infinite level of truncation where 
all the higher derivatives are taken into account. In table 1 we present results obtained 
through the procedure of successive truncations and through our numerical solution for 

/o\ 

N — 3. We give the values of k, A, u\ (n) for the scaling solution and the critical 
exponents p, v. We observe how the results stabilize as more p-derivatives of Uk(p) at 
p = k and the anomalous dimension are taken into account. The last line gives the results 
of our numerical solution of eq. (|5|). By comparing with the previous line we conclude 
that the inclusion of all the p-derivatives higher than u^\k) and the term ~ p in the 
“threshold” function of eq. @ generates an improvement of less than 1 % for the results. 
This is a lot smaller than the error induced by the omission of the higher derivative 
terms in the average action, which typically generates an uncertainty of the order of the 
anomalous dimension. In table 2 we compare our values for the critical exponents with 
more accurate results obtained with other methods (such as the e-expansion, summed 
perturbation theory at fixed dimension, lattice calculations and the 1/iV-expansion). As 
expected p is rather poorly determined since it is the quantity most seriously affected by 
the omission of the higher derivative terms in the average action. The exponent v is in 
agreement with the known results at the 1-5 % level, with a discrepancy roughly equal to 
the value of p for various N. 

In conclusion, the shape of the average potential is under good quantitative control for 
every scale k. This permits a quantitative understanding of the most important properties 
of the system at every length scale. Our investigation is not restricted to the behaviour 
near the phase transition on which we have concentrated here because it is the most 
difficult to handle numerically. Also the initial form of the potential does not have to be 
of the quartic form of eq. (|T9|) . Arbitrary general short-distance potential can be studied. 
For example, the tricritical point for the transition to a first order phase transition can be 
investigated with our numerical methods. Furthermore, the flow equation is well defined 
for arbitrary continuous d. This will permit an explicit check of the validity of the e- 
expansion for more detailed quantities characterizing the equation of state. 
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Tables 



ft* 

A* 

( 3 ) 

Ui 

V 

V 

a 

6.57 x 10" 2 

11.5 



0.745 

b 

8.01 x 10" 2 

7.27 

52.8 


0.794 

c 

7.86 x 10" 2 

6.64 

42.0 

3.6 x 10" 2 

0.760 

d 

7.75 x 10" 2 

6.94 

43.5 

3.8 x 10” 2 

0.753 

e 

7.71 x 10" 2 

7.03 

43.4 

3.8 x 10" 2 

0.752 

f 

7.64 x 10" 2 

7.07 

44.2 

3.8 x 10" 2 

0.747 


Table 1: The minimum k of the potential Uk(p), the derivatives A = u"(n ), uf\n) for the 
scaling solution, and the critical exponents rj and u, in various approximations: (a)-(e) 
from ref. || and (f) from the present letter. N = 3. 

a) Truncation where only the evolution of k and A is considered and higher derivatives of 
the potential and the anomalous dimension are neglected. 

/o\ 

b) k, A, u k (k) are included. 

/o\ 

c) k, A, u k (k) are included and rj is approximated by eq. fllUp. 

d) with five parameters: k, A, u^\k), u^\k) and rj. 

e) as in d) and in addition u^\k), u^\k) are estimated. 

f) The partial differential equation (||) for u' k (p) is solved numerically and p is approxi¬ 
mated by eq. rn- 


ii 















N 

V 

V 

1 

0.643 

0.6300(15)“ 

0.6310(15) fe 

0.6305(15)“ 

0.044 

0.032(3)“ 

0.0375(25) b 

2 

0.697 

0.6695(20)“ 

0.671(5) b 

0.672(7)“ 

0.042 

0.033(4)“ 

0.040(3) 6 

3 

0.747 

0.705(3)“ 

0.710(7) b 

0.715(20)“ 

0.038 

0.033(4)“ 

0.040(3) 6 

4 

0.787 


0.034 


10 

0.904 

0.877 d 

0.019 

0.025 d 

100 

0.990 

0.989 d 

0.002 

0.003 d 


Tabic 2: Critical exponents v and ij for various values of N. For comparison we list 
results obtained with other methods as summarized in |T0| and 0: 

a) From summed perturbation series in fixed dimension 3 at six-loop order. 

b) From the e-expansion at order e 5 . 

c) From lattice calculations. 

d) From the l/TV-expansion at order 1/N 2 . 
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Figures 

Fig. 1 The evolution of u' k (p) as k is lowered from A to zero. The initial conditions (bare 
couplings) have been chosen such that the scaling solution is approached before the 
system evolves towards the symmetric phase with u' k { 0) > 0. N = 1. 

Fig. 2 The evolution of k, A and rj for the solution of fig. 1. 
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